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Abstract 

Quantum mechanical phases arising from a periodicahy varying Hamihonian are considered. 
These phases are derived from the eigenvalues of a stationary, "dressed" Hamiltonian that is able 
to treat internal atomic or molecular structure in addition to the time variation. In the limit of 
an adiabatic time variation, the usual Berry phase is recovered. For more rapid variation, non- 
adiabatic corrections to the Berry phase are recovered in perturbation theory, and their explicit 
dependence on internal structure emerges. Simple demonstrations of this formalism are given, to 
particles containing interacting spins, and to molecules in electric fields. 
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I. INTRODUCTION 



Any quantum mechanical system in a stationary state accumulates a dynamical phase 
over time proportional to the energy of that state. To determine energy differences, based on 
phase differences accumulated between two such states, is the basis of Ramsey spectroscopy; 
the workhorse of high-precision measurement. For this reason, small effects that can add 
spurious phase shifts must be understood and kept under control. 

One such variety of spurious phases arises when the states of interest are not strictly 
stationary. If the total Hamiltonian H{t) has an explicit time dependence, then this depen- 
dence will generate an additional phase evolution. For example, precision spectroscopy of 
trapped ions must contend with the fact that the ions are in motion, and experience varying 
ambient fields during the course of their orbit. Likewise, ultracold atomic samples confined 
in an optical dipole trap experience, in principle, an oscillating electric field due to the laser 
field that provides the trapping potential. 

In the case of a Hamiltonian with a slow, periodic time dependence H{t+T) = H{t), Berry 
|l| has given a famous description of the additional phase. Berry's original treatment requires 
that the period r be far larger than any other relevant time scale of the system, and thus finds 
an "adiabatic" phase shift. This shift is largely independent of the detailed way in which 
the Hamiltonian varies with time, and leads to an elegant geometric description of the phase 
. Extensions to this formalism have considered the next-order corrections 
if the rate of change of the Hamiltonian is not strictly adiabatic 
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A more general Floquet theory has also been advanced, which allows one to consider the 
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16| . In addition, the ideas have been 



17|, gauge structure [l8|, the Quantum Hall 



effect of overtones of the fundamental period r 
extended to particles with dynamic properties 
Effect [3], and to relativistic effects using the Dirac equation [20|. 

Thus far, applications of the Berry phase have mostly considered the effect of the time- 
dependence on quantum mechanical particles without internal structure, although atoms 



with two or several levels have been considered [16|, l21l . l22l |. However, the job of precision 
spectroscopy is precisely to reveal this internal structure. Corrections to Berry's phase aris- 
ing from degrees of freedom internal to an atom or molecule is our concern in this article. 
To establish a concrete formalism for this, we will consider a particular case, namely, a dia- 
magnetic or electrically polar species in the presence of a magnetic or electric field, whose 
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FIG. 1: The axis of rotation with laboratory-fixed coordinates f?, C} ^ weh as the field coor- 
dinates defined by {x,y,z}. The field direction rotates about the C-axis with angular frequency 

direction precesses on a cone with an angular frequency uJr (see Fig. [I]). The system evolves 
in time according to the field variation, combined with whatever intrinsic Hamiltonian gov- 
erns the particle's internal structure. The internal structure dictates regimes of linear and 
quadratic Zeeman (Stark) shifts with respect to the applied magnetic (electric) field. 

A main point in deriving the non-dynamic phase in this situation is to recognize the 
periodicity of the driving field. By analogy to the periodic driving of a near-resonant laser 
field applied to a two-level atom, we consider "field-dressed" states of the Hamiltonian in 



the spirit of quantum optics 
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231]. This viewpoint effectively counts the energy of the 
atom itself, plus that of the photons of frequency Uj. arising from the driving field. The 
additional energy shift due to the rotating field is then equivalent to the ac Stark effect in 
optics. By constructing the complete Hamiltonian in this way, we are able to accommodate 
the particle's internal structure. We are also able to consider arbitrary rates of rotation, 
not just those that are adiabatic with respect to the particle's Hamiltonian. Nevertheless, 
in what follows we will focus primarily on results for low rotation rates, to better draw 
analogies with the usual adiabatic phase. 

This article is structured as follows. In Sec. II we work out the general transformation 
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from a Hamiltonian with an explicit time-rotating field, into an effective dressed Hamiltonian 
whose eigen-energies yield the shifted energies. We illustrate this transformation first with 
a simple two-level atom, then generalize it to an arbitrary atom or molecule. In Sec. Ill 
we briefly re-visit a structureless particle with arbitrary total spin, showing that our results 
reduce to Berry's in the limit of slow rotation. Sec. IV illustrates the application of the 
method to a particle composed of two interacting spin-1/2 objects, to show most clearly the 
effect of their fine structure on the Berry phase. Finally, Sec. V considers a simple case of a 
dipolar molecule in a rotating electric field, to assess the influence of molecular end-over-end 
rotation on the phase. 

II. DERIVATION OF THE DRESSED STATES 

This section lays the ground work for the formalism. It is composed of two subsections. 
The first works out the derivation for a spin— 1/2 system in a rotating magnetic field, and will 
illustrate simply and cleanly the basic idea. The second subsection contains a generalization 
of this same derivation to explicitly account for an arbitrary quantum mechanical object 
with total spin j. 

A. Spin— 1/2 interaction 

The most elementary of quantum mechanical objects is the spin— 1/2 particle. Its inter- 
action with a time- dependent magnetic field is 

H{t) = -fl:^-B{t). (1) 

There are two coordinate systems that we will find useful, the lab-fixed Cartesian coordinates 
labeled and a coordinate system {x,y,z) rotating with the magnetic field, where 

z = B (see Fig.[T]). In the lab-fixed frame the field rotation is given by its components 

B{t) = {B^{t),B,{t),Bc{t)) 

= B {sm9rCOSLJrt,sm9rSmLJrt,cos9r) , (2) 

where (j) = uJrt. The field makes an angle Or with respect to the lab-fixed C, axis, which is 
also the axis it rotates about. In other words we identify a vector ujr = oJrCi fhe axis defining 
the field's rotation. 
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The spin-field interaction has its own characteristic frequency, 

- Hm- B = — - — a ■ B = -ujlo- ■ B, (3) 

in terms of the Pauh matrices a and the (m-independent) Larmor frequency ui = q^bB. 
g is the usual (7-factor; e.g. g ~ 2.0023 for the simple case of the spin-1/2 electron. The 
eigen-energies of the non-rotating spin-field interaction are ±h ul/^^,. By causing the field 
to rotate, we expect to find apparent energies that are different from these. From this point 
forward we shall work in atomic units and set /i = 1. 

In this simple example, we can easily cast the Hamiltonian in matrix form, using the basis 
of spin functions = ±1/2), with quantization along the axis. These states themselves 
are then independent of time, which would not be the case if we used the eigenstates referred 
to the (moving) magnetic field. Using the explicit forms for the Pauli matrices, we get 

^^^^ ^ujlI cos(^,) sin(^,) e--^* \ 
2 I sin(^^) e*^"-* - cos(^,.) j 

While this basis seems cumbersome it is quite useful for the dressed-state calculation we 
wish to perform. 

We begin in the most general way by writing an ansatz wave function in this basis: 

\m) = . . h (5) 



me 



where we leave the frequencies uj± unspecified for now. We note a very similar development 
has recently been carried out by the authors of Ref. [24] for a two-level system. The constants 
will be chosen later to cancel any time dependence in the system. The wave function 
satisfies the time-dependent Schrodinger equation (TDSE) 

t^ = Hitm. (6) 

Inserting the ansatz ([5]) into the TDSE yields the equation of motion for a and f3. This new 
equation introduces terms which amount to a kind of effective Hamiltonian due to rotation. 
Moving these to the RHS of the TDSE, we get 



+ ^ COS 9r-uj+ ^ sin ^^,e-^K+'^— '^+)* 
\ I I ^sine^e*(^'-+^-^'^+)* -!^cos^^-cj_ 



(7) 



Formally, this is still a time-dependent Hamiltonian unless we use the freedom in choosing 
the frequencies u± to get rid of this dependence. We could achieve this by setting 

UJr-\- OJ- — OJ-^- (8) 

There are many ways to accomplish this, but a particularly symmetric and appealing one is 
to let 

^± - ±y , (9) 

or, even more to the point, 

^m<; = ixir (10) 

for each angular momentum projection along the rotation axis. This last statement is 
very general, and will motivate our choice of wave function in the spin-j case below. 

In any event, we are presented with a formally time-independent Hamiltonian, dressed 
by the rotation: 

-f^dressed = -f^nr + -f^r, (H) 

where the magnetic-field part, in the laboratory-fixed |m^) basis, looks like 




2 I sin(6'r) — cos 

and is exactly the same as having the magnetic field tilted at an angle 9r with respect to 
the quantization axis in the direction of the x-axis, as specified by = ujrt — 0. Rather, the 
vestiges of rotation show up in the effective term 

/ 
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= — diag(m^ ujr). (13) 



v 



-1 
+1 

Now that we have a time-independent Hamiltonian we write 

where A is the effective, or dressed, eigen-frequency of the "stationary" state defined by 
i?nr + Hy-. A is, in other words, the eigenvalue of i^dressed^ 

A± = ±^-y/a;£ - 2uL0Jr cos{9r) -\- oSl- (15) 



exp(-iAt), (14) 



Thus the apparent energy difference, measured by the phase difference accumulated during 
a hold time r, would be (A+ — A_)r = — rather than the (1^+1/2 — '^-i/2)t = — V^- 
that would measure the intrinsic energy splitting. The difference between these two we 
attribute to a Berry-like phase in each state defined as (where we assume the dynamical 
phase is unaffected by the rotation) 

l± = <P±-^±: (16) 

where (f) [ip) is the total (dynamic) phase of the system. This removes the dynamical phase 
if that would have been accumulated in the absence of rotation. Hence 7/r represents the 
error introduced by the field's rotation into a measurement of the state's energy. 

The dressed eigen-energies are actually independent of whether or not the time r refers 
to one period of the rotating field. To connect explicitly to Berry's phase, let r = 2-71 /ur 
be one period of the field's rotation. Berry's phase results when the field rotation is slow, 
i.e. when the frequency of the rotating field is small compared to the Larmor precession 
frequency of the spin. If we assert that Ur <^ uji, then the spin is expected to precess around 
the instantaneous field direction, and to follow this direction during the rotation of the field. 
In this limit the dressed-state eigenvalue of one state, say A+, is 

A+-f -yCos(^.). (17) 
In a period r the additional phase accumulated by this state is (see (fT6ll ) 

^ n -^{2n)cos{9r). (18) 

Here we have used the freedom to add vr to the phase, which amounts to an unobservable 
overall change of sign of the dressed eigenstate (oq Po)- 

The ordinary dynamical phase = ujlT is canceled out in flTH]) leaving only the terms 
due to the slow rotation of the field. This remainder can be written 

7(0) = ^(27r)(l - cos(0,)) = mAQ, (19) 

where 

An= d<p I sin(^)c/^ = 27r(l - cos(0^)) (20) 
Jo Jo 
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is the solid angle subtended by the vector describing the direction of the rotating field. In 
other words, this additional phase 7^°-* is exactly Berry's phase mAn, for a particle with 
spin projection m = 1/2 along the (^-axis. Similarly, for the m = —1/2 state there arises an 
additional phase —{1/2)AQ. 

We have cast the Hamiltonian H in (ITT]) in the basis where angular momentum projection 
along the rotation axis Ur = C is a good quantum number. Yet, in the limit of slow rotation 
that is particularly important, the relevant states are those where m is quantized along the 
field axis. We should therefore cast H in the basis of states \sm) quantized along the field 
axis, which we call the i-axis (see Fig. 1). 

Using the ansatz our dressed Hamiltonian represents a Hamiltonian if^r that is 
already diagonal in this basis, and reads 
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1 
-1 



(21) 



To cast the rotational piece in this basis, we make explicit the notation |s m^) for spin 
states quantized along the rotation axis C; and |s m) for spin states quantized along the field 
axis z. The rotation from to z is accomplished through the Euler angles R = (0,6*^,0). 
The transformation from |s rrii^) to |s m) is given by 25 1 



\s m) = D{R)\s m^) 



(22) 



(23) 



The Wigner rotation matrices have a simple explicit form: 



(24) 



where 



,1/2 _ J , _ 
^1/2,1/2 ~ '^-1/2,-1/2 — ^ Y 



jl/2 



"-1/2,1/2 



''^1/2-1/2 — sm I 2 



(25) 
(26) 
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Using the fact that is diagonal in the rotation basis with values —rriQUOri it is easily verified 
that, in the rotating-frame basis, 



if. = ^ ..o,.., ^^^^ 




Therefore, the Hamiltonian matrix referred to the field axis is 

H = H^, + H, (28) 
^-^^cos^. -^^sin^, 

(29) 



2 — 

2 — -r- 2+2 



-i^sin^. -^^^ + ^cose,. 



And this matrix gives has the same eigenvalues as f[T^ . as is expected when one merely 
performs a unitary transformation on the system. 

In this basis, the Hamiltonian is already diagonal in the absence of rotation. Thus for a 
small rotation it is nearly diagonal, and the eigenvalues are easily estimated in perturbation 
theory. Indeed, Berry's energy follows immediately from the diagonal perturbations in this 
matrix: it is —{uJr/'2) cos{6r) for m = +1/2, and +{uJr/2) cos{9r) for m = —1/2. 

To summarize, for a Hamiltonian of the form H{t) that we have been dealing with, the 
rotation-dressed energies are given by the eigenvalues of the time-independent operator 

-f^dressed = H{t = 0) - tU^-Sf = H{t = 0) - Ur {cOs{dr)Sz - Sm{6r)Sx) ■ (30) 

For more general systems incorporating internal structure, all that is required is to include 
the appropriate structure in H{t = 0), as we will now see. 



B. A general spin-j system 

Having provided the groundwork by working out the simple spin- 1/2 system, we now 
proceed with a general derivation for the spin-j system. A system such as this can be 
described by the following Hamiltonian 

H{t) = Ho-fl-f{t), (31) 

where J-'{t) is an external field, electric or magnetic, that acts on an appropriate moment 
of the atom or molecule. JF rotates on a cone at frequency Ur and tilt angle Or just as B did 
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in the previous section and depicted in Fig. HI Here Ho is a Hamiltonian in the absence of the 
apphed rotating field. It can be used to describe the hyperfine elements of an atom or it can 
be a detailed molecular Hamiltonian that includes such items as rotation, spin-spin, nuclear 
spin, or lambda doubling. The Hamiltonian can equally be represented in a basis referred to 
the axis of rotation or to the instantaneous field axis. Later we will take the instantaneous 
field axis, just as we did above. We will, as in the previous section, begin by quantizing 
along the axis of rotation. We take this structured object and place it in a rotating field JF; 
this can be either electric or magnetic provided there is an electric or magnetic dipole that 
interacts with the field in the usual way, i.e. it is a scalar interaction of two vectors. 
To work with this Hamiltonian, it is again convenient to pick two basis sets: 

|(/€)jmf) 

\{K)jm,) (32) 

Because j is the total of all relevant angular momenta, its projection onto an axis is unam- 
biguously defined as in the lab fame and rrij in the rotating frame, as above. Here k is 
a shorthand notation for all the other quantum numbers required to specify the state. 

To deal with the explicit time dependence of the field rotation, we will expand into the 
lab basis first, and will make the same trial wave function that was motivated above; 

\m)= E C^«'.',-^e-^'"c-*|(«:')jH)- (33) 

We have explicitly included a time dependent phase factor with phase rrit^uJr, which is akin 
to the spin-1/2 case, cf. ( JTOjl . Taking the time derivative for the TDSE and projecting onto 
a particular state, gives 

{(K)jm^\i^ = {tC^,j,m, + m^cu.a,,, e-^™c-'*. (34) 

As for the internal Hamiltonian Hq, it may or may not be diagonal in our basis, but it does 
not depend on any external field. Therefore it can be represented in a basis where it is 
diagonal in m^, whereby 
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To treat the field interaction, we use the language of tensor algebra, and express the 
spherical components of T in the lab frame as an explicit rotation from T in the rotating 
frame (whose z axis is, of course, defined by the instantaneous direction of T itself): 

L 

= J^Vl;iuJrt,9r,0). (36) 

JF is the magnitude of the field, and q is its spherical projection in the rotating frame. But JF 
defines this frame so only the values of g = will contribute. D is a Wigner rotation matrix. 
In a similar manner, the dipole moment fl is determined by its spherical components such 
that 

= -^$^(-i)>.Pi:oM,^r,o) 

= -^^/i,e-^-^*c?i,o(^,). (37) 

I — I 

This uses the explicit expression for T> in terms of a little-d function [25]. 

Just as we treated the internal degrees of freedom in Hq we must now treat the field 
interaction. 



{{K)jm(;\ - fl ■ J^\ip) = 2^ (-l)'((«:)j"^cl/^J('«')j"K) ^ 

C.',',m'diA9r)e'-^'^-'-"'<^. (38) 

Piecing together the different parts, and multiplying through by e*'"^'^''*, we arrive at the 
TDSE for the coefficients C: 

iCK,j,m^ + i^rmt^C^j^rn^ = ( (/t) jm,^ | | (/t') j'm^) C^/ j/^^^ - 

K',j' 

^ X] (-i)''((/t)j"^cl/^J(«;')i'K)*^'^'.i'.K ^ 

k' j'm'^L 

dl^o{9r)e''^'-'^-'-"'c+"'^\ (39) 

— * 

Now, we have not specified what field JF is, nor which structural degrees of freedom are 
involved in making the dipole fl, and it does not matter. All that matters is that fl is 
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a vector, in which case the Wigner-Eckhart theorem apphes 25|. In the total angular 
momentum basis, we must have 

((z/)jmd/xJ(z/OjH)oc I ^ ^ ^ I , (40) 

where the proportionality constant involves the reduced matrix element. Then the con- 
servation of angular momentum implies that m'^ — = —l. However, this immediately 
removes the time-dependence in the exponential term in (l39l) . In fact, as alluded to earlier, 
this statement says that any angular momentum imparted by the rotating field must be 
accounted for in the projection m,^. 

In some cases it will prove more useful to keep track of the individual spin components 
rrii^i separately. For instance, suppose there were two angular momenta, m^i and m^2'- we 
would have two projection terms that would each evolve as e*'"ci'^''* and e^'^f^'^r*^ Terms in 
the Hamiltonian which describe the interaction of the two spins are of the form si ■ S2 for 
which we find scales as 



{si m^i S2 m^2|si ■ S2\s[ m'^^ ^ m'^^) 



oc 




The proportionality involves a reduced matrix element. By the conservation of angular 
momentum we find that m^i — m'^^^ = —(777,^2 — ^(2) ^"^^ phase factor is still canceled out. 
In fact, for any such interaction between two spins, the conservation of angular momentum 
forces the time dependence to cancel out. 

With the time-dependence removed, Eq. (15^ reduces to the Schrodinger equation for a 
non-rotating field tilted at an angle 6r from the rotation axis, just as in the spin-1/2 example 
in the prior subsection. This introduces an additional term on the left of fl39l) . which is moved 
to the RHS and interpreted as an effective Hamiltonian. Thus if Hq is presented in the basis 
rrij) diagonal with respect to the field, then the matrix to be diagonalized is 

-f^dressed = Hq — fl ■ T — UrTTLQ 

where 777^ is the eigenvalue of the operator. As before, we can rotate this Hamiltonian 
from the Tn.^ basis to the rrij basis. Since Hq does not depend on the either 777^ or 777^-, it is 
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unaffected by this rotation. In the frame of the instantaneous field where mj is the good 
quantum number, we can write the dressed Hamiltonian as 



-f^dressed — Hq — [1 ■ J-' — LJr ■ j , 

= Hq- fl-f - uJr{cos{9r)jz " sin(6'^)j^). 



(43) 



This dressed Hamiltonian is the main result of this article. This Hamiltonian has been 
previously formulated in NMR studies 



261]. In the following sections we will apply it to 



several elementary cases of interest. 



III. PURE SPIN-s SYSTEM 



As the simplest application of the general method beyond the spin-1/2 particle, we con- 
sider in this section a structureless particle of arbitrary spin s, as was considered in the 
original formulation of Berry This spin interacts with a magnetic field that rotates at 
an angle 6r with respect to the axis of rotation. Using the result from (H3|) . this system is 
described by the Hamiltonian 



-f^dressed = ^LSz " UJr{cOs{6r) S ^ - Sm{6r)Sx), 



(44) 



where Ul = Osf^B^ is the m-independent Larmor precession frequency and Qs is the (^-factor 
for the spin-s. For this section, we have reverted to the usual notation s and m for the spin 
and it's projection onto the instantaneous field axis. 

For this structureless particle, the Hamiltonian is represented by a (2s + 1) x (2s + 1) 
tridiagonal matrix, in the basis of states \sm). This matrix is explicitly given by 



H 



m a 6(s, m) 

b{s,m) {m — 1) a 6(s,m — 1) : 

b{s,m — l) '•■ '•■ 6(s, — m + 1) 

: 6(s, — m + 1) — (m — 1) a b{s,—m) 

b{s, —m) —m a. 



(45) 



where a = lul — cos(6r) and 6(s,m) = {l/2)^/s{s + 1) — m{m — 1) ujrsm{6r). Appendix 
A sketches a derivation of the eigenvalues of this matrix, which are 



Am = TTlJ Uj\ + Uj"^ — 2uJr-UJL COs(6'r), 



(46) 
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where m takes on the values — s, . . . , +s in integer steps. 

The usual Berry phase is obtained in the adiabatic limit where cu,. <^ cj^, in which case 



Am ~ muji — mujr cos(0,.). (47) 

The magnetic field completes one rotation in a time r = 27r/uJr- In this time the spin 
accumulates a dynamical phase = muJiT. Beyond this, it acquires an additional phase 
7m = 0m — V^m, whcrc 0m = Am^, givcu to lowcst ordcr by 

= -m27r cos(^^) ^ 27rm(l - cos{e,)). (48) 

In the final step, we use the fact that adding 2 tt m (where m is either integer or half-integer) 
amounts to adding an unobservable phase of ±1 to the system. We see that the phase 7m'' 
accumulated is exactly that given by the result of Berry; m times the solid angle subtended 
by the rotation. We can also extend this solution to regimes of non-adiabaticity. The first 
order correction in uir/uj^ is 

7« =27rm^sin2(^,), (49) 



which has already been identified elsewhere jsl, 14]. Based on our explicit formula, we can 



extract corrections to any desired order, at least for fields undergoing the simple motion 
in Fig. [B This additional phase can be expanded to any desired order in the adiabatic 
parameter Ur/uJi. For example, the second and third-order 7^'^'-' corrections are 



7^2) = 27rm-^cos(^,)sin^(^,) (50) 



■'L 



= 27rm-^(3 + 5cos(2e,))sin2(0,) (51) 



,3 

Using the general dressed formalism, the limit of fast field rotation can also be described. 
The phase 7^ (after one field period) can be approximated in this limit {ur ^ ul)'- 

7m ~ 27rm — 27rm— cos(6'r), (52) 

and the first term is unobservable. In this case the dominant energy, as manifested in the 
phase, is the photon energy due to the time-periodic field. On top of this, the magnetic field 
interaction itself makes a small correction. This is clearly not the appropriate limit in which 
to perform precision spectroscopy, since small uncertainties in the field rotation rate would 
dominate the observable Larmor frequency. 

14 



The structureless spin problem can be solved analytically for arbitrary rotation rates since 
the spin precesses about an effective magnetic field in the rotating coordinate system j2^. 
The effective magnetic field makes an angle with respect to the positive ^-axis, 6**, satisfying 



1 - 2^ cos(^,.) + 



, 2 



cos(^J - ^ 

cos(^:) = (53) 
l-2^cos(^„) + (^) 

This angle smoothly transitions from 9* ~ 6r for Ur <^ ujl to 9* ^ n for Ur ^ ujl- Equiv- 
alently, the angular deviation, A9r = 9* — 9r, between the effective magnetic field and the 
true magnetic field can be described through 



sin(A6'r) 



2^ cos(^,) +{^'^ 



1 - ^ COsi9r) 

cos{Mr) = (54) 
y^l-2^cos(^.)+(^)' 

where this deviation smoothly transitions from A9* ~ for a;^ -C cui;, to A9* ^ n — 9r for 
The total phase shift, 0^ = V^m + 7m5 accumulated during one revolution of the magnetic 



field can 
where 



3e broken up into a dynamic contribution, and a geometric contribution, 7^, 



26 



28 



29| 



if^ = 27cm— cos{A9r), 

= 27rm(l-cos(e:)). (55) 

ip*^ and 7^ are slightly different than the phases defined in Eq. ( fT6l) . In this definition, both 
the dynamical phase and the geometric phase acquire non-adiabatic corrections. In the limit 
of very slow rotation cu^, we have that — > Grouping terms together and simplifying 
provides the total accumulated phase 



2vrm I l + i^Jl-2^cos(^,)+ ( ^ 




27rm l + Wl-2— cos(^^)+ — . (56) 
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The total phase shift arrived at in this geometrical way is the same as the dressed state 
derivation — [271 / ujr)\mi where is defined in (1461) — apart from a factor of 27rm, which 
is unobservable for integer or half-integer values of m. This is a semi-classical geometric 
procedure that yields the non-adiabatic corrections to Berry's result. 

For the above calculations, we have taken pi = —gsf^sS and Ur to have a positive sense 
about the (^-axis. For the single revolution phase shifts, 7m\ terms proportional to odd 
powers of Ur/oJi change sign when the g-factor changes sign, while terms proportional to even 
powers of uj^/uol change sign when the sense of rotation changes sign. For the expressions 
defining the angle, 6**, and angular deviation, A^,., changing the sign of the g-factor (sense 
of rotation) directly changes the sign of ojl {oJr)- 

A spectroscopic measurement would involve finding the energy difference between two 
states with different values of m, with difference Am. In a Ramsey-type experiment, this 
measurement seeks to measure the phase difference uJiAmT = A(f. In a rotating field, how- 
ever, the experiment will produce a measurement of A0 = A■m^yuJ1 + uj"^ — 2ujLUJr cos(6'r)r, 
and thus will introduce an error. This error is given by the difference A7 = 7m — 7m' = 
A0 — Ay9, and is plotted in Fig. IIIII as a function of rotation rate. The different curves 
represent different values of the tilt angle Or- 

IV. STRUCTURED SPIN-J SYSTEM 

More generally, atoms and molecules are composite objects made of individual spins, 
which are moreover coupled together to create fine or hyperfine structure. For example, 
alkali atoms couple the electronic and nuclear spins into a total hyperfine state. The resulting 
angular momentum structure will have a bearing on the non-adiabatic corrections to the 
geometric phase accumulated. 

As a simple illustration of our formalism, we consider a composite particle composed of 
two spin-1/2 objects. This example goes beyond the structureless particle often envisioned 
by the usual Berry theory. The dressed Hamiltonian is given by 

^^dressed = ^ijlz + ^232z + Aji ■ J2 " UJr{cOs{0r) J ^ - sin(6',.) J^r) , (57) 

— * — * 

where uJi = QiHb^ is the Larmor precession frequency of spin jf, J is the vector sum of ji 
and j2] and A is parameter that governs the splitting between levels J = and J = 1. 
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10-1 10° 
a)r/ajL 



10^ 10^ 



FIG. 2: The extra phase accumulated due to the rotation of the field. In the limit of very fast 
rotation, uJr ^ ujl, the system accumulates a phase of 47r, which is unobservable. In this limit, the 
states are best represented by projections onto the axis of rotation. The various lines represent 
values of Or between 7r/2 (bottom line) and (top line) in steps of vr/lO. As can be seen, when 9r 
is zero, there is no measurable phase shift, since there is no enclosed solid angle. 



The Hamiltonian fl571) can be represented by a 4 x 4 matrix, in the basis { | (ji 72) J Mj) } 

m |)oo),i(H)ii), l(H)io)' l(H)i -1)} 



/ _3A 
' 4 



dressed 



v 





UJl — UJ2] 





V2 





— Ur C0S( 

sm(9r) 



\{UJI - UJ2) 



V2 



A 

4 



V2 





sin( 



\ 







(58) 



^ sin(^^) ^-uz + ujr cos(^^) / 



where uz = + (^2); is the average of the individual Larmor frequencies. The first 

item to note is that if the two spins have identical Larmor frequencies, uji = 002, then 
this Hamiltonian is equivalent to that of a spin-0 particle and a spin-1 particle that are 
independent of each other, there is no coupling between the two states. Each would then 
evolve according to the previous section on pure spins. This would be the case for the singlet 
and triplet excited states of the helium atom, (ls2s)^''^5' state, for example. However, should 
these spins be different from one another (such that uJi 7^ 002) then coupling corrections arise. 

The ordinary adiabaticity criterion specifies that the rotational frequency Ur be small 
as compared to the Larmor precession frequency ujl, which in this example is given by 
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uji. However, now it becomes also necessary to specify whether the Larmor frequency itself 
is large or small compared to the splitting A between adjacent J- levels. This is because 
the Berry phase arises from a correction to the eigenvalues of the dressed Hamiltonian 
relative to the non-rotating Hamiltonian. It is therefore worthwhile to cast the non-rotating 
Hamiltonian in the basis in which it is as diagonal as possible. In the following sections 
we treat the two limits separately. Since our emphasis here is on the Berry-phase limit, 
we consider only the limit where lOr <^ ci^z, where the rotation rate of the field is small 
compared to the Larmor frequency. The resulting phase shifts are of course implicit in the 
theory, however. 

A. Weak magnetic field, lo^ <C <C A 

In the low-field limit, but assuming that each Larmor frequency lOi is still far larger than 
the rotational frequency cUr, we can write down expressions for the energy quite simply. 
Note that in the absence of rotation, the leading-order energy shift is the sum of the Larmor 
frequencies themselves, i.e., (1/2) (cui + co'2)Mj = oj^Mj. Leaving this correction on the 
diagonal to break the degeneracy of the J = 1 level, we now treat as perturbations the 
difference (1/2) (cji — 002) and the rotation rate uo^. 

Doing so, the leading-order correction due to rotation of the field is given by the diagonal 
terms in (l58ll that contain the rotation rate Ur- This correction is the usual Berry phase 
found above, 

-i^^\\J Mj)) = -27rMjcos^„ (59) 

and it depends on the atomic state only through the total projection of angular momentum 
Mj. Thus the ordinary Berry phase in the limit of zero rotation rate is still intact, and is 
independent of the internal structure. 

However, the higher order corrections do depend on this structure. To leading order in 
the rotation frequency Ur, we find a correction to the Berry phase in the |(| |)1 1) state: 

V'.(|n))^2.|;-..^(«(i.,^^). (60) 

The first term is the usual first order correction for a structureless particle (cf. (H9l) ). with 
the replacement of u;^ by uz- This should be expected since the energy splitting between 
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the two states is given by cjz, and thus is what must be overcome by the rotating field that 
couples together the differing projections. The second term in the parentheses depends on 
how strongly the rotating field couples states with differing total J, as manifested by uj^—oj^- 
This new correction arises from 4*^-order mixing in perturbation theory, it is nevertheless 
linear in the adiabaticity parameter uo^/uoz- A similar expression is found for the |1 — 1) 

As is expected from the pure spin case, this state picks up an overall negative sign. However, 
due to slight changes introduced by the structure, we find a slightly different correction to 
the second term in parentheses. In fact, we can write down an expression that encapsulates 
the first order (in ojr/oJz) correction as 

M.>, ^ 2.M.^ si„^(«.) (l + M. ^il^^ . (62) 



While the first term is exactly of the form in (H9i) . the second term describes how the 
distant \ J = 0,Mj = 0) state affects the accumulated first order phase 7mj; namely that 
the quadratic Zeeman shift in the two Mj = levels distorts the system such that the 
I J = 1, Mj = 1(— 1)) state is affected more (less) by the | J = 1, Mj = 0) state. 

The two states with Mj = do not acquire a geometric phase at lowest order in Ur, which 
is appropriate. In this case the leading order perturbation to the dressed Hamiltonian is 
= ±[(l/2)(co'i — uj2)]'^/A, which denotes the quadratic Zeeman shift already present in 
the non-rotating system, and which does not contribute to the Berry phase 7. The quadratic 
shift after a period t = 27t /ur is the dynamical phase the Mj = would nominally acquire. 
To the first order in Ur in which there is a correction to the | J = 0, Mj = 0) state arises in 
4*'^-order perturbation theory. It is given by 

7(^)(|00)) = 2vrg^sin^(g.) , (63) 



2A ' A{i-{-^y 



where as always the superscript "1" denotes a correction linear in Ur- Here we find a term 
that appears similar to the Mj = 1 states, with the exception that it occurs in an Mj = 
state. We have introduced a new energy scale into the problem by adding A and this allows 
the 1 0) state to acquire a first order Berry phase. However, the strength of this phase is 
reduced by a term proportional to the ratio of the quadratic Zeeman shift in the lower level 
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to the spin-spin energy splitting. Given our assumptions, this term — while linear in ujr — 
is a product of multiple small parameters, and is generally smaller than 7*^°^ for Mj = ±1 
states. 

For the case of the | J = 1, Mj = 0) state, there is also a 4*'^-order correction, but it takes 
a very different form. We find, after much algebra, 

7«(|10)) = -27r-^sin2(^.)^±^. (64) 

ZUJz l^Z 

This is very different from the |0 0) state correction in (jHS])- The important energy scale 
is the linear Zeeman shift ujz- Eqs. (1631) and ( |64l) carry an important insight; the energy 
scale responsible for higher-order Berry phases is different for the two Mj = states. The 
dominant scale in the |0 0) is the spin-spin splitting A. In the |10), the dominant energy 
scale is the linear Zeeman shift. To first order in ujr, there is a correction to the Mj = 
states that, while similar to the shift in the \Mj\ = 1 states, is reduced in magnitude. This 
reduction is due to the structure, the structure that provides a quadratic Zeeman shift in the 
Mj = states. For the lower (upper) level, the correction depends on the relative strength 
of the quadratic Zeeman shift to the spin-spin splitting (linear Zeeman shift). In the regime 
considered, both of these contributions are very small. The same ideas apply to the F = 
and F = 1 hyperfine states of Hydrogen, where the magnetic field is coupling states of the 
same parity. Briefly, 7(1) is influenced by "nearby" Mj = ±1 states for the |10) level, and 
comparitively less influenced by the "far away" Mj = ±1 levels in the |0 0) state. 

It is instructive to examine these results for different cases of individual Larmor frequen- 
cies. In the case where both particles experience the same Larmor frequency in a field, 
uJi = UJ2, then these first-order corrections reduce to the usual first-order corrections for a 
structureless spin-1 particle, as in Eqn. fH9|) . and the additional Mj = pieces are zero as 
well. In another limit where one Larmor frequency dominates the other, say ui ^ co'2, then 
7^^) for the \Mj\ = 1 states reduce to the first order correction of the dominant spin alone, 
reflecting the fact that the weaker spin is coupled to the stronger one and gets dragged along 
for the ride. This happens, for example, in the F = 1 hyperfine ground state of the hydrogen 
atom, where the nuclear g-factor is far smaller than the electron g-factor. 



20 



B. Strong magnetic field, cji , a;2 S> A S> 



In the other hmit, where the magnetic field is large compared to the splitting between 
adjacent J-levels, it is more useful to construct the dressed Hamiltonian in an alternative 
basis. Namely, the non-rotating Hamiltonian is more nearly diagonal in the independent-spin 
basis \ ji mi j2 1^2) , where the four Zeeman energies -E'mi,m2 are given simply by miU!i + m2U!2'- 

E_l_l = ~{UJI+LJ2) 

E_i^i = ~{uji-uj2), (65) 

as appropriate to this Paschen-Back limit of the Zeeman effect. This is an example of 
the aside in Sec. II where we make the ansatz (a, P) e.~'^'^^<^rt^-im2uirt _ fpjjg remaining 
Hamiltonian, which includes the rotation of the field and the spin-spin interaction, is recast 
as follows 

/ £^1 1 + f - uJrCos{er) f sin(^^) sin(^^) \ 



dressed 



fsin(e,) ^i,-i-f fMOr) -t 



^ sin(^^) E_i__i + f + a;,. cos(^^) -f sin(^^) 

V -fsin(e,) -f -fsin(e.) ^-i|-fy 

(66) 

Once again, we can immediately read the Berry-phase contribution from the diagonal com- 
ponents, as 

7^°Hlii "^1 32 m2)) = -27r(mi + 1112) cos{9r) 

= -27rMjcos(^^). (67) 

The phase accumulates due to the individual spins separately, as expected when the spins 
interact weakly with each other compared to their interaction with the field. 

This independent accumulation of phase leads to a different interpretation of the Mj — 
states: in the limit of small magnetic field compared to the spin-spin energy splitting, we 
had attributed this to an Mj — projection while here, we can attribute this to mi = ±| 
accumulating ±{u!r/2) cos{9r) extra energy and the m2 = t| accumulating =p(a;r/2) cos{9r)- 
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To cement this idea even further, we find there are two independent first-order contributions 
to the first-order non-adiabatic correction 7*^^^ computed first neglecting A: 

mi j2 7712)) = 27T mi-^sin^ (Or) + 2n 7712-^ sin^ {6r) (68) 
= l['\\jimi))+^i'\\j2m2)). (69) 

This is exactly the contribution one would expect from two independent spins following a 
rotating field. When the spins are anti- aligned, or mi = —m2, we find no correction at this 
order. It is worth noting that this perturbative expansion breaks down if Ur ^ uji. Thus, 
should the rotation rate be fast with respect to one of the Larmor frequencies, but not the 
other, then the measured phase difference cannot be treated perturbatively in this regime. 

The explicit effect of internal structure, manifested in the splitting A, appears as a next- 
order correction: 

,a.^2.™.i^.n^(«„)(i.(A) 

Of course, there are many routes by which 4*'^-order perturbation theory can affect this state. 
We only give one route which produces a phase shift proportional to Ur after one period of 
oscillation. There is a structure correction for each non-adiabatic spin that depends on the 
relative strength of the spin-spin splitting to the other Larmor frequency. Thus, we require 
that the spin-spin splitting be small compared to each of the Larmor frequencies in order 
to make this expansion. Again, this is a more restrictive condition on adiabaticity than is 
usually employed for two independent spins. 




V. POLAR MOLECULES IN A ROTATING ELECTRIC FIELD 



Molecules bring yet another degree of freedom to the picture, namely end-over-end rota- 
tion with eigenstates M^v). In addition, if the molecule is polar, it has an electric dipole 
moment that can be acted upon by a rotating electric field. In this section we will consider 
only diatomic molecules, and only one of a fairly simple structure, to illustrate how our 
formalism applies to them. The lowest-order Berry phase was worked out recently in this 
system 22], but the higher-order corrections are implicit there as well. 
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For the sake of illustration we choose the simplest of diatomic molecules, a molecule 
with no hyperfine structure. In a rotating electric field this system is described by a Hamil- 
tonian of the form 

^^dressed = BN"^ - flm " £ ~ UJr{cOs{9r) - Sm{9r)N^), (72) 

where N is the end-over-end rotational angular momentum of the molecule, /im is the electric 
dipole moment of the molecule and S is the electric field strength. Since we are in the frame 
of the electric field and the electric dipole moment points along the molecular axis, we do 
not have any couplings of or A, where A is the projection of total angular momentum 
onto the internuclear axis. For S-molecules, this means there are no couplings to excited 
electronic states by the applied electric field at the low fields we consider. 

For simplicity, we consider here the coupling only between the = and A^ = 1 rotational 
levels of the molecule, assuming weak coupling of rotational states due to the electric field, 
i.e., fim^ ^ B. The formalism can of course be extended to arbitrarily large A^ values as 
needed. It is nice to note that this formalism has an atomic analog: the ^5*0 and ^Pi states 
of noble gas and alkaline-earth atoms have opposite-parity and are coupled by the Stark 
interaction. This approach gives the corrections for states of opposite parity coupled by 
the Stark interaction. The dressed Hamiltonian reads, in the basis {|A^Mjv)} = {|00), |1 — 

i),|io),|i + i)} 

/ -j^f^mS \ 



dressed 



2B — Ur COs{9r) "^^r Sm{6r) 



^fJ'rnS -^ujrsm{9r) 2B -^w^ sin( 



(73) 



\ -^ujr sin(6'r) 2B + Ur cos(6'r) / 

Note the the electric Hamiltonian is off-diagonal in the basis of parity eigen-states. In the 
absence of the perturbation uj^, this Hamiltonian appears to have a complete degeneracy 
among the three states with A^ = 1. In the magnetic field case above, this degeneracy was 
broken by the linear Zeeman effect acting on the diagonal matrix elements. To achieve the 
same feat here, we must account for the off-diagonal mixing due to electric field. Note the 



similarity of this procedure to that of Vutha and DeMille 22 1 . 



We first diagonalize the Mat = subspace, using the mixing angle 6 defined by 

tan(5) = = -X, (74) 
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with the usual eigenvectors (cos(5/2), sin((5/2)) and (— sin((5/2), cos(5/2)). The exphcit val- 
ues in terms of the parameter x are 



.S\ VT+?+l 



In terms of this mixing angle the transformed Hamiltonian, with electric-field-dependent 
terms on the diagonal only, reads 

I B{l-VlTx^) -^sin(f)sin(^,) sin(f ) sin(ft,) \ 

-^sin(|)sin(^^) 25 - cj^ cos(0^) ^ cos(f ) sin(^^) 



dressed 



^ cos(|) sin(ft,) B{l + VlTx^) ^cos(|)sin(^,) 

\ — ^sin(|)sin(e^) ^ cos (| sin (ft,) 25 + w,. cos(^^) / 



(77) 



In the limit that a; <^ 1 we see that the diagonal terms for the two = states are merely 
the quadratic Stark shift = ±{nmSy/QB. From here we can read off the ordinary 

Berry phase from the diagonal perturbations linear in Ur, yielding the usual 

7(°)(l^^iv)) = -27rM;vcos(^,), (78) 



where by is meant the appropriate eigenstate of the field-mixed Mn = states 22] . This 
diagonalization removed the degeneracy of the M^r = level with the Mjy = ±1 levels of 
the = 1 subspace. However, the degeneracy among the Mn = ±1 states still exists. 

Having quasi-broken the degeneracy in the = 1 levels, we can now evaluate the first- 
order adiabatic correction term using standard second-order perturbation theory. It is ev- 
ident that both states |A^, M^r = ±1) experience the same additional phase at this order 
(due to their degeneracy in the non-perturbing Hamiltonian), given by 

7«(|1, ±1) = -TT u;rsm\9r), (79) 



y(/i„£:)2/ 

where an expansion in the small parameter x has been applied. Requiring this to be a small 
correction identifies the adiabaticity criterion for this situation. If this case were analogous 
to the magnetic field case, we would only be concerned about the magnitude of uJr with 
respect to fim^, which is the stand-in for the Larmor frequency in this case. However, 
Eqn. (1791) suggests a slightly different criterion, namely UrB <ti {iim^Y must hold in order 
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to recover the simple leading-order Berry phase. To understand the origin of this criterion 
we look at the term in the large parentheses in fl79|) . It is the inverse of the Stark energy in 
the absence of field rotation for the upper level. We can rewrite (17^ as 

^^'\\NMm = ±1)) = -27r^sin2(^,), (80) 

and we recover a form reminiscent of the pure spin case (cf. (H9!) ). where the Larmor 
frequency, ul, is replaced by Eg, the quadratic Stark shift. In order to be an adiabatic 
correction, it is immediately evident why UrB <^ {jj,m£Y must hold; the rotation rate must 
be small compared to the energy splitting in that level. In this case, the splitting is quadratic 
in electric field and therefore a secondary energy scale — the rotational level splitting or 
internal structure — must come into play. 

By similar reasoning, we can arrive at the 7^^-' corrections for the Mjy = states as well. 
That of the lower level is given by 

7«(|iV~0,0)) = -27r||sin2(^^„)^, (81) 

where we find the requirement that uj^ ^ must hold. In this case, our assumptions 
clearly support this adiabatic criterion since we are in the regime Ur <^ fJ^m^ <^ B. This 
correction is linear in Ur, but suppressed by the ratio of the Stark energy to the rotational 
constant of the molecule. Physically, this is because this state is far removed from the 
"degeneracy" in the = 1 levels. It is the electric analogue of the weak magnetic field limit 
of the spin-spin interaction. The correction for the upper level is given by 

7(^)(|iV~l,0)) = 4vr^sin2(^.). (82) 

This correction is in the opposite direction to and twice that of the |A^Mjv = ±1) states 
because the ~ 1,0) state is influenced by the two states, M^v = ±1, that are below 
it in energy. This is in contrast to the coupled spins in a magnetic field case, where each 
Mm = ±1 contributed equally in magnitude but opposite in sign. This is due to the lack of 
any linear Stark shift in the Af/v = ±1 levels. 

It is evident that polar molecules in a rotating electric field are quite similar to magnetic 
dipoles in rotating magnetic fields. There is an energy splitting in comparison with which 
the rotation of the field must be small to ensure adiabaticity. If there is a shift in energy 
that is linear with the applied field, then the rotation rate must be small compared to this 
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energy. However, if the energy scales quadratically with the apphed field, the rate of rotation 
must be small in comparison to the energy shift in the field. Thus, the internal structure 
is quite important in regimes of quadratic field shifts and introduces different adiabaticity 
requirements on ujj. in terms of the applied field and internal structure. 



VI. CONCLUSIONS 



The dressed-state formalism is a natural way to treat quantum mechanical objects subject 
to time-periodic driving such as we have considered here. It allows for the inclusion of 
arbitrary internal structure of the object, and still reveals the exact dressed eigen-energies 
at arbitrary rotation rates. In the limit of slow rotations, it also reduces, as it must, to 
the usual geometric Berry phase. Because it includes the structure of the atom or molecule 
considered, however, it is also able to shed light on the infiuence of this structure on non- 
adiabatic corrections to the geometric phase. It is therefore expected to be a powerful tool to 
)e used when analyzing high-precision spectroscopic data in the presence of periodic driving 
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The treatment herein has considered only the simplest case of a magnetic or electric field 
whose direction precesses uniformly about a given axis. It is to be expected, however, that 
this treatment is yet more general, and that dressed states for arbitrary periodic driving 
Hamiltonians could be constructed, at least numerically. It could, for example, be combined 
with the Floquet analysis that has been used previously for unstructured particles 
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APPENDIX A 

In Sec. Illlt we had a Hamiltonian of the form 

H = ujlSz — uJr{cos{6r)Sz — sm{6r)Sx)i (Al) 
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which looks hke 



H 



ai Ci 

ci 02 C2 : 

C2 C„_2 

• C„_2 On-1 C„_i 



(A2) 



y . . . c„_i an- j 
In our system, not only is this matrix symmetric, it contains the following added symmetry: 
a\ — —ttn, 02 = — On-i, etc. Also, the coupling coefficients q follow a similar property: 
Ci = c„_i, etc. These properties are key to simplifying the eigenvalues of the tridiagonal 
matrix in our case. 

Eigenvalues of an n x n tridiagonal matrix are given by the roots of the polynomial Pn, 
defined recursively by 

Po(A) = 1, 

Pi (A) = (oi-A), 

P2(A) = (02 - A)pi(A) - CiPo(A) 



Pn{X) = (a„ - A)p„_i(A) - 4_iP„_2(A). 



In our problem, the constants are defined as 



dm = m{u;L - cos(6'r)), 

S 1 S 
-m q m — q 



(A3) 

(A4) 
(A5) 



Thus, the symmetry pops right out. 

A simple example is the case of 5" = 1/2. Here we find that the characteristic polynomial 

is 

P2(A) = (A - ^a)(A + ia) - Qfe^ = (A^ - \{a' + b')) (A6) 

where a — {cul — u!rC0s{9r)) and b — u!rSm{9r). For the case of 5" = 1 we find a similar 
equation (after simplification) 



(A7) 



P3(A) = X{{X-a){X + a)-2i^—b] ) 
= A(A^ - {a' + b')) 
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For integer values there is always a diagonal element that is 0. As is evident, this has the 
same form as flA6|) with the added piece of A multiplying everything yielding an eigenvalue 
of 0. In addition, Eq. (1A7P is scaled by a factor of 4 from from Eq. (lA6p . thus making the 
the eigenvalues a factor of 2 larger. This is because the value of m in flATD is twice as large 
as the value of m in (lA6p . 

We can generalize the characteristic polynomial to a very simple expression due to the 
added symmetries. It is given by 

P2S+i= n (A2(i-^™-in>o)-m2(a^ + 6')). (A8) 

m=(mmin>0) 

The Kronecker (5-function in (lASP is to insure that in the event mmin = there is only one 
eigenvalue A = 0. 
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